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1. INTRODUCTION 


The purpose of this paper is to present an analytic model of an imaging sensor system so that: (1) sen- 
sor performance predictions can be made; (2) design tradeoffs and sensitivity analyses can be rapidly 
performed; and (3) insight into various aspects of imaging sensor performance can be obtained. The 
model is applicable to image sensors which operate from the visible through the thermal infrared spec- 
tral regions. 

The design of sensors for remote observation of the Earth from polar orbiting satellites takes about a 
decade and occurs in five different stages. These stages are: 

• Definition of Scientific Requirements. During this stage, scientific working groups formulate 
the scientific requirements (e.g., spatial, temporal, and spectral resolution; measurement ac- 
curacy, etc.) 

• Preliminary Design. During this stage a “paper design” is developed. (Determination is made 
of such parameters as f-number, detector size, and number of detectors) that permit the sen- 
sor to meet the scientific requirements (e.g., signal-to-noise ratio, noise-equivalent delta tem- 
perature, or noise-equivalent delta reflectance). 

• Feasibility Studies (Phase B). During this stage, engineering feasibility is established without 
regard to optimization of the design. 

• Design Studies (Phase B). During this phase, the design is optimized, and an in-depth analy- 
sis is performed on each subsystem (e.g., optics, focal plane, cooler, electronics, mechanical 
systems) and a credible cost estimate is produced. 

• Flight Hardware Design, Development, Test and Integration into the Space Platform 
(Phase CjD). During this stage, flight hardware is designed, developed, and tested to prove 
that it meets the specification, and is integrated into the space platform. 

Only Preliminary Design (the second stage) is addressed in this document. Derivations or references 
are given for all the equations to make it easy to change the theory as required in future applications. 

The spectral range is limited to 0.4 to 15.0 pm which is generally appropriate for studies of the Earth 
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and its environment. The types of scanners include the “pushbroom” and two different kinds of 
“whiskbrooms.” A substantial portion of the analysis presented herein has been incorporated into a 
self-documented Lotus 1-2-3 spreadsheet. 

Analytic (as opposed to statistical) methods are used in the model. Analyses are carried out at very 
low spatial (and therefore temporal) frequencies in order to simplify the computations. High spatial 
frequencies are used only to determine the MTF of the sensor. Carrying out the analyses at low spat- 
ial frequencies also avoids the necessity of working in the frequency domain which generally involves 
fourier transformations or convolutions in the time domain, and would therefore make the sensor 
model very complex. To further simplify the analyses, it is also assumed that the sensors have narrow 
spectral bandwidths so that those parameters that are spectrally dependent (e.g., detector responsivity) 
can be reduced to a constant. 

The model assumes that the sensor is a linear system from the optical input through the electronic 
signal processing. This assumption is satisfied if (1) the optical system is not dominated by diffraction 
effects; (2) incoherent detection methods are employed; (3) the various noise sources are additive in an 
rss sense; (4) the imaging process is spatially invariant; and (5) the electronic signal processing consti- 
tute linear operations. These assumptions are all valid for the types of imaging sensors that the model 
is presently being applied to. 

Appropriate scene radiance levels must be assigned in order to assess sensor system performance. In 
the model described in this document, tables are included which allow a user to determine radiances 
at the top of the Earth-atmosphere system as a function of wavelength. These tables apply to the 
visible and near ir spectral regions. The origin of the tables is discussed in Appendix B. For the 
thermal infrared region, radiances are directly computed. 

Section 2 addresses the power incident on the detector during an observation interval and the signal 
coming out the detector. These are written as functions of the irradiance at the detector and as a 
function of the scene radiance. 
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Section 3 provides a description of the model used for the scene radiance. 


Section 4 derives the equations for detector dwell time for four different scanning configurations which 
include: a spinning 45-degree scan mirror, a spinning paddle mirror, a rocking scan mirror, and a 

Unear array-pushbroom. 

Section 5 is devoted to the major noise sources associated with visible and infrared detectors. 

Section 6 presents definitions and derivations of various “figures of merit including Noise Equivalent 
Delta Reflectance (NEAp) and Noise Equivalent Delta Temperature (NEAT). 

Section 7 develops the Modulation Transfer Function (MTF) for the optical aperture, the detector aper- 
ture, satellite motion, charge diffusion, charge transfer, and satellite jitter. 

The Appendixes treat many of the equations in a tutorial manner. 
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2. POWER AND SIGNAL 


When viewing the scene, the power (flux) into the detector is given by 


$ = 


E(A) dX 


[W], 


where Ap = 
E(X) = 
X 

AX = 


the area of the detector (cm 2 ); 

scene spectral irradiance at the detector (W/cm 2 -pm); 
center wavelength of spectral band (/im); and 
X 2 - Xj = spectral bandpass (;um). 


( 2 - 1 ) 


The signal out of the detector is given by 

/• X-j 


S = A 


D 


E(A)R(A) dX 


[A], 


where R(A) is the detector current responsivity and is given by 


( 2 - 2 ) 


where q 

V 

h 

c 


q?? 

R(X) = X [A/W], (2-3) 

h c 

= the charge of an electron = 1.60 X 10“ 19 [coul] ; 

= the detector quantum efficiency [nd] (assumed to be constant over the spectral 
bandpass AX); 

= Planck’s constant = 6.63 X 10~ 34 [Wsec 2 ]; and 
= the speed of light = 3.00 X 10 10 [cm/sec]. 


Substituting Equation (2-3) into Equation (2-2) gives 


S = 


A D qr? 

he 


rh 

I E(A)A dX 


[A], 


J X, 


(2-4) 
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Equation (2-4) may be written in terms of the spectral photon irradiance as 

rh 


s = A d qr? 


E'(A) dX 


tA], 


JX 


(2-5) 


where the scene spectral photon irradiance E'(X) at the detector is given by 


A o 

E'(X) = E(X) [p/sec-cm -jum] 


he 


( 2 - 6 ) 


The signal out of the detector is given by 


*1 

S' = S [e] 

q 


(2-7) 


where the integration time t[ is given by 

t t = Ct D [sec], (2-8) 

where C is an input constant and t D is the sensor dwell time, which is described in detail in Sec- 
tion 4. 


Substituting Equation (2-5) into Equation (2-7) gives 

/•Xo 


S — t| Aj-j T\ 


E'(X) dX [e] . 


J Xi 


(2-9) 


The scene spectral irradiance at the detector is related to the scene spectral radiance by the following two 
equations given in terms of watts and photons, respectively. (See Appendix A.) The first equation is 

r 7T 

E(A) = L(A) [W/cm 2 -pm], (2-10) 

4f 2 N 

where t q = the optical transmittance [nd] from the sensor aperture to the detector, and 
f^ = f-number of the optics [nd] . 

The second equation is 


E' (A) - 


r 0 7r 


4r 


L'(A) 


[p/sec-cm 2 -pm] 


N 


(2-H) 
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Also, for completeness, note that 

X i 

L'(X) = L(X) [p/cm-sec-sr-pm] . (2-12) 

he 

For the visible and shortwave infrared (SWIR) wavelengths, it is assumed that the spectral bandpass 
AX is small and that the scene spectral photon irradiance E'(X) varies slowly over AX. In this case, 
the integration in Equations (2-1) and (2-4) can then be approximated by E'AX, where E' is the 
average value of E(X) over the spectral bandpass AX. 
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3. SPECTRAL RADIANCE 


In this section, we will show how to obtain the spectral radiance for the visible and SWIR spectral 
wavelengths (0.4 to 2.2 pm) and the infrared spectral wavelengths (2.2 to 15 pm). 

3.1 Visible and Shortwave Infrared (SWIR) Spectral Radiance 

A table of values of the scene spectral radiance at the sensor aperture is given in Appendix B as a 
function of p, 6 Z , and X 

where p = Earth’s surface reflectance (nd), and 
6 Z = solar zenith angle (deg). 

The scene spectral radiance L(X) is obtained from the tables by trilinear interpolation at the de- 
sired values of p, d z , and X. Although the values of scene spectral radiance in the tables were 
computed for a nadir viewing sensor, we assume that the values are independent of viewing angle. 
(See Equation C-4). 


3.2 Infrared Spectral Radiance 

The total scene spectral radiance in the infrared is given by 

L(X) = t a B(X,T s ) + e A B(X,T a ) [W/cm 2 -sr-pm] 
where r A is the atmospheric transmittance for an optical depth 5, given by 


t A 


e 


-5 


(3-2) 


(3-3) 


and where e A is the atmospheric emissivity given by 

•a - • - (W) 

The quantities B(X,Tg) and B(X,T A ) are the spectral radiances of the Earths surface at temperature Tg 
and the atmosphere at temperature T A , respectively, and are given by Planck’s equation (Hudson, 

1969, p. 35) for a blackbody at temperature T, 
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B(A,T) 




[w/cm 2 -sr-/im] 


(3-5) 


where Hudson’s equation (2-8) has been divided by n to convert to radiance, and where 
C, = 2?rhc 2 = 3.74 X 10 4 [W-qm 4 /cm 2 ] ; 

h = Planck’s constant = 6.63 X 10 -34 [W-sec 2 ]; 

c = speed of light = 2.998 X 10 !0 [cm/sec]; 

C 2 = hc/k B = 1.44 X 10 4 [jum-K]; 

k B = Boltzmann’s constant = 1.38 X 10 -23 [W-sec/K] ; and 
T = Blackbody temperature [K]. 


The total scene spectral photon radiance is given by 

L'(A) = r A B'(A,T S ) + e A B'(A,T A ) [p/sec-cm 2 -sr-/im] , (3-6) 

where the surface spectral photon radiance B'(A,T<-) and the atmospheric spectral photon radiance 
B'(A,T a ) are found by evaluating the following equation (Hudson, 1969, p. 38) at and T A , 
respectively: 

C 'l 1 

B'(A,T) = — [p/sec-cm 2 -sr-/im] (3-7) 

“• (-©-0 

where 

C', = 2ttc = 1.88 X 10 23 [p-sec _I -cm' 2 -gm 3 ] 
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4. DWELL TIME 


As the satellite moves in orbit, it images along scan lines perpendicular to the ground track. Let 
the extreme ends of the scan lines on the Earth subtend an angle 0 at the satellite. This angle 
is called the field-of-view (FOV). See Figure 1 for an illustration of the geometry involved and refer 
to Appendix C for a discussion of the relations between the parameters shown in the figure. The 
number of angular resolution elements n F is given by 

n F =— [nd] (4-D 

n a 

and (as can be seen from Figure 1) the instantaneous angular field of view a (geometric only) 
is given by 

a = ^l [rad ] (4-2) 

where d s = detector width (mm), and 

f = focal length of the optical system (mm). 

A scan line can be formed by a spinning mirror, a rocking mirror, or a linear array. Each of these 
are discussed in the following sections. 

4.1 Spinning Mirror 

Two types of spinning mirrors are discussed in this section: the 45-degree scan (faceted) mirror and 

the paddle scan mirror (Figures 2a and 2b). The axis of rotation for a 45-degree scan mirror is paral- 
lel to the sensor’s optical axis, and a change in the angle of rotation of 9 will produce an equal 
change in the line-of-sight angle. The axis of rotation for a paddle scan mirror is perpendicular to 
the sensor’s optical axis, and therefore, a change in the angle of rotation of 9 will result in a 26 
change in the line-of-sight angle. 

The dwell time for a spinning mirror is given by 

t D = — [sec] , (4-3) 

n E n f 


9 



d s = Width of square detector (/rm) 

f = Focal length (cm) 

H = Satellite height (km) 

R e = Radius of Earth (km) 

R = Earth-to-satellite distance (km) 
s 


S d = Maximum scan angle distance (km) 
a = Instantaneous angular FOV (deg) 

0 = Swath width FOV (deg) 

= Ground swath angle (deg) 

S w = Swath width (km) 


Figure 1. FOV and IFOV Geometry 
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SIDE VIEW 


FRONT VIEW 



’NOTE: AS MIRROR ROTATES THROUGH AN ANGLE 6 THE LINE OF SIGHT ALSO ROTATES 
THROUGH AN ANGLE d 

Figure 2a. Scan Mirror Geometry for 45-Degree Faceted Scan Mirror 


PADDLE SCAN FLAT MIRROR 



’NOTE: AS THE MIRROR ROTATES THROUGH AN ANGLE d THE LINE OF SIGHT ROTATES 
THROUGH AN ANGLE 20 

Figure 2b. Scan Mirror Geometry for Paddle Scan Flat Mirror 
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where rip = the number of facets, 

! > I for a 45-degree scan mirror, and 
= 1 for a paddle scan mirror. 

The active scan time t^, which is that part of the scan mirror period during which data are acquired, 
is given by 

t A = «t M [sec], (4-4) 

where the scan mirror period t M is the time for the spinning scan mirror to make a complete revolu- 
tion, and the scan efficiency k is given by 


K 


n ( -Q 

2 ^ 


[nd], 


(4-5) 


where 4 = 1 for a 45-degree scan mirror, 

f m{ 

(= 2 for a paddle scan mirror. 

During each scan, an area on the ground is imaged. The extent imaged along the ground track at 
nadir (Figures 3 and 4) is 

d c = n f n D aH [km], (4-6) 

where H = the satellite height [km], and 

nj)= the number of detectors per spectral band. 

Also, during each scan, the satellite moves a distance d m measured along the ground track, which is 
given by 


d m { M V SUB 


where V SUB is the speed of the subsatellite point along the ground track and is given by 

(]i)V 2 


V SUB ~ R e 


(R p + H) 3 / 2 


where 


li = GM = 3.98603 X 10 5 (km 3 -sec' 2 ) 

where G is the universal gravitational constant, and M is the mass of the Earth. 

The overlap factor B is given by 
d c 

B = — — [nd], 

d m 


(4-7) 


(4-8) 


(4-9) 
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SATELLITE 



d m = Ground distance moved per scan period (m) 

d c = Extent imaged along ground track per scan period (m) 

0 = Swath width FOV (deg) 

H = Satellite height (km) 

S = Swath width 

w 

V SUB = Subsatellite velocity (km/sec) 

Figure 3. Ground Track and Scanning Geometry 
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FACET #1 


j, 

FACET #2 


FACET #3 


I 

FACET #4 

J 


DETECTOR #1 








DETECTOR #2 








DETECTOR #3 








DETECTOR #1 








DETECTOR #2 








DETECTOR #3 








DETECTOR #1 



■ 





DETECTOR #2 








DETECTOR #3 








DETECTOR #1 








DETECTOR #2 








DETECTOR #3 



■ 







n f = 4 = Number of facets on scan mirror 
n D = 3 = Number of detectors per spectral band 
H = Satellite height (km) 

S w = Swath width (km) 

a = Instantaneous angular FOV (as in Fig. 1) (deg) 
d c = Extent imaged along ground track per scan (m) 

Figure 4. Example of Scanning Geometry for a 45-Degree Faceted Scan Mirror 
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where B = 1 for contiguous coverage, B > 1 for overlap, and B < 1 for incomplete coverage (under- 
lap). In terms of the percentage of overlap S 0 we may write 


B = 1 + 


100 


[nd] 


(4-10) 


Using Equations (4-7), (4-9), and (4-6), one obtains 

n t -n D aH 


l M ~ 


BY 


sec 


SUB 


However, from Equations (4-1), (4-3), (4-4), and (4-5), tp may be written as 


l D ~ 


a 


2ttF 


M 


t M [sec] . 


(4-11) 


(4-12) 


From Equation (4-12) it can be seen that t D is independent of the field of view 0, and therefore, 
changing the swath width does not change the dwell time. Substituting Equation (4-11) into (4-12) 
yields 

*i f ii D c 


D 2it bv sub f m 


[sec] . 


(4-13) 


4.2 Rocking Mirror 

During a complete cycle, the rocking mirror rotates back and forth through an angle that covers the 
FOV 0 twice— once in each direction. For this case, the dwell time is 

t n = — — — [sec] (4-14) 

n S n E 

where n<j = 1, for imaging in the forward scan direction only, and 

n<j = 2, for imaging in both directions. 

The active scan time t A is given by Equation (4-4). However, in this case, the scan efficiency k de- 
pends on the sensor design parameters; e.g., mirror turnaround time, value of n s , and mirror inertia. 
The scan period T M in this case is the time for the rocking mirror to do a forward scan and retrace 
to its initial position. 
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During each scan period an area on the ground is imaged. The extent imaged along the ground track 
at nadir (Figure 3) is 


dc = n s n D a! ^ [km] • 

Use of Equations (4-7), (4-9), and (4-15) results in 




n s n D «H 


BV 


[sec] 


SUB 


(4-15) 


(4-16) 


Using Equations (4-4), (4-14), and (4-16) gives 

[sec] . 




k n D aH 


n E BV SUB 


(4-17) 


4.3 Linear Array 

In this case, no scan mirror is used, and the number of angular resolution elements n E is equal to the 
number of detectors in the cross-track direction. The dwell time for a linear array is, therefore, 


«H 

t D = — [sec], (4-18) 

V SUB 


and the field of view © is given by 

0 = n E ex [rad] . (4-19) 
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5. NOISE 


The total sensor noise is composed of detector noise and system noise. The detector noise is com- 
posed of a number of different noises that depend on the material composition of the detector and 
whether it is a single detector or is used in an array. 


The total sensor noise is given by 

n tot = ^ n ~det + n2 sys 

where N DET = root-mean-square detector noise, and 
N S ys = root-mean-square system noise. 

The following sections present a discussion of these various noise sources. 



(5-1) 


5.1 Visible and SWIR Detector Noise 

Here we will describe the various types of noise encountered in visible and SWIR detectors. 

5.1.1 Visible and SWIR Detector Noise Sources 

5. 1.1.1 Photon Noise 

Photon noise or shot noise is due to the random arrival of photons at the detector. Because the inci- 
dent photon flux follows a Poisson distribution, the photon noise is given by (Levi, 1968, p. 153) 

N p = (S') 1 / 2 [e] , (5-2) 

where the signal S' is the mean number of electrons produced by the photons arriving at the detector 
and is given by Equation (2-9). 

5. 1.1. 2 Bulk Trap Noise 

Bulk trap noise occurs in CCD focal plane arrays and arises from the random trapping and emission 
from interface or bulk states (Dereniak, 1984, p. 243) and is given by 

N BT = ( M g k B TA D n SS 8n 2 ^ 1/2 (el, (5-3) 

where M = the number of gate transfers; 

o 

n ss = the density of surface states; 

A d = the detector area; 
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k B = Boltzmann’s constant; and 
T = Temperature [ K ] . 

Typical values of N BT are: 

• Surface Channel CCD— 1000 electrons 

• Buried Channel CCD— 100 electrons. 

5. 1.1. 3 Output Amplifier Noise 

The output amplifier noise is associated with the amplifier that buffers the signal from the focal plane 
and is generally a metal oxide semiconductor field effect transistor (MOSFET) of a given transconduc- 
tance. An expression that can be used to compute this noise is (Dereniak, 1984, p. 244) 


N 


8C 0 2 k B TAf \ 1/2 


OA 


3< l 2 gm 


[e], 


(5-4) 


where C. 


the output capacitance [pfarad]; 
the electrical bandwidth [Hz] ; 
the transconductance of the MOSFET [mhos]; 
the charge of an electron [coul] ; and 
Boltzmann’s constant. 

5. 1.1.4 Dark Current Noise 

Dark current or thermal generation noise is associated with charge carriers that are thermally generated 
to bring the CCD potential well into thermal equilibrium. Dark current root-mean-square (rms) noise is 
given by (Honeywell, 1986, p. 5-26) 


o 

Af 

q 

k B 


N 


2j dc a d 4 i 


1/2 


DC 


[e] 


(5-5) 


where = dark current density at temperature T [a/cm 2 ]; 


V D 


= detector area [cm 2 ]; 

= integration time [sec]; and 

= electron charge [coul]. 


The dark current density gi yen by 
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(5-6a) 


J DC = a M T3 ex P 



where E 


g 


“M 


silicon band gap = 1.12 eV; 

Boltzmann’s constant = 8.62 X 10” 5 [eV/K]; 

temperature [K] ; 

a material-dependent carrier recombination factor; for silicon r? M = 2; and 
a material-dependent factor, typically = 1.1 X 10 6 [A/cm 3 -K 3 ] 


5. 1.1. 5 Johnson (Thermal) Noise 

The thermal motion of electrons in a resistor gives rise to voltage fluctuations across the resistor leads. 
These fluctuations are known as Johnson or thermal noise. The noise current is given by (Dereniak, 


1984, p. 39) 

) l/2 

,A], 

where Af = the effective bandwidth of the circuit [Hz] and 

R = the resistance [12], 



(5-7) 


It follows that the noise in electrons is given by 



5. 1.1. 6 Schottky Noise 

Electrons in the semiconductor of an M-S (metal-semiconductor) junction may overcome the potential 
barrier to reach the metal and produce a noise current. The noise current is called Schottky barrier 
noise and is given by (Yang, 1978, p. 130) 


A d r c t 



[A], 


where R c = 47rqmk 2 /h 3 = Richardson’s constant = 120 [A/cm 2 -k 2 ]; 


(5-9a) 
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/V D 

k B 

T 

q 


detector area [cm] ; 

work function = 0.0354 [eV], for PdSi:Si diodes; 
Boltzmann’s constant = 8.62 X 10 -5 [eV/K] ; 
temperature [K] ; and 
electron charge [coul]. 


The Schottky noise current can be converted to electrons and is given by 


N DCS “ 



[e]. 


(5-9b) 


5. 1.1. 7 Charge Transfer Noise 

Charge transfer or transfer inefficiency noise is associated with CCD structures and occurs because of 
the random amount of charge lost by a signal upon transfer and the amount of charge introduced to 
a signal upon entering a well. The noise N^j associated with a single well (Dereniak, 1984, p. 242) 
is given by 



where e = the transfer efficiency [nd] - 


(5-10) 


If the number of detectors is n' D , and the number of phases to transfer the charge is n p , then the 
total number of wells is n p n' D . Hence the total charge transfer noise is given by 


N 


CT 


= ^en^npS 


1/2 


[e] 


(5-11) 


5.1.2 Visible and SWIR System Noise 
5. 1.2.1 Quantization Noise 

The quantization noise Nq is given by (Montgomery, 1978, p. B-l) 

S SAT 


Q 12 1 / 2 2 ^ 


[e], 


(5-12) 


where Q is the number of bits used in the analog-to-digital (A/D) converter. From Equation (2-9) 
one obtains for the visible and SWIR 
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(5-13) 


S SAT _ t I A D T?E SAT t e l> 

where S'^j is the signal that would result if the detector were receiving the saturation irradiance. 
This is the flux that produces a signal level that just causes the A/D converter to saturate. 


The saturation irradiance E'g A -p is given by 

E' sat = E M S F [p/sec-cm 2 -^m] , 


(5-14) 


where, from Equation (2-11), E' 


M 


e 'm = 


T 0 n 


L' 


4P 


M 


N 


is given by 
[p/sec-cm 2 -jum] . 


(5-15) 


When the saturation factor Sp is multiplied by the maximum expected scene irradiance E' M at the 
detector, an irradiance E'^j will be produced that will just saturate the A/D converter. 


5. 1.2. 2 Other System Noise 

When system noises are from unknown sources, they are designated as other system noise. 


5.2 Infrared Detector Noise 

5.2.1 Infrared Detector Noise Sources 

The detector noise is composed of two parts, photon noise and other noise. It is given by 

/ \ 1/2 

N DET = ^ N P 2 + N OD 2 j ’ (5-16) 

We will calculate the photon noise Np on the basis of a cold shielded and cold filtered detector. The 
other detector noise Nqj->, will be estimated from laboratory values of D*. We will estimate the other 
detector noise from D* first. 

5. 2. 1.1 Other Detector Noise 

Consider the detector that will be used in the sensor (same area and electrical bandwidth), but with 
background temperature and viewing angles identical to those used in the laboratory measurement of D*. 

The other detector noise Nq D is given by 
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n od =(nee 2 -n pl 2 ) 


1/2 


[e] . 


(5-17) 


The other detector noise Nqq is assumed to be independent of the level of cold shielding and cold filter- 
ing. The photon noise under laboratory conditions is 

v l/2 


N PL ~ ( l l F L ^ 7 A D j" B ^ T BGL ! ^ d ^ 


[e], 


(5-18) 


where is the detector cutoff wavelength, T BGL is the background temperature in the laboratory, 
B'(T BGL > IS defined by Equation (3-7), and F L is the view factor in the laboratory (F L = 1 when 
viewing 180°, or 2 tt steradians). In general 


F L = sin2 (-y) [nd] 
where is the full-cone view angle. 

The detector noise equivalent power is given by 

(A d AO 1 / 2 


(5-19) 


NEP - 


D* 


[W], 


(5-20) 


where A D = 
D* = 


detector are [cm 2 ]; and 

the laboratory value of specific detectivity [cm-Hz^/W] . 


The number of noise equivalent electrons NEE is given by 

Rt t 

NEE = NEP [e] (5-21) 

q 

and the effective noise bandwidth is given by 

Af = — [kHz], (5-22) 

2tj 

where |3 is the ratio of noise bandwidth to information bandwidth, and tj is the integration time given 
by Equation (2-8). 
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5. 2. 1.2 Photon Noise (Infrared) 

The photon noise is given by (Levi, 1968, p. 153) 


N p = (r') 1/2 [e], (5-23) 

where the total number of electrons produced by the scene and background is given by 


II 

D 

-3 

- r^2 

E'(X) dX + 

(e'bg 00 

[e], 

(5-24) 


. Jx, 

Jx', J 




where X] to X 2 is the spectral bandpass of the scene photon spectral irradiance E'(X), which is given 
by Equation (2-11), and X', to X' 2 is the spectral bandpass of the background photon spectral irrad- 
iance E' bg (X), which is governed by the cold filter. 


The background spectral irradiance E' BG (X) is given by 


e'bg <*) = 


S2 e Q + bG - ^ T 


CF 


B'(X,T bg ) [p/sec-cm 2 -pm] (5-25) 


and the emissivity of the optics is given by 


e 0 = 1 - r 0 fndl , (5-26) 

where B'(X,T BG ) = the background spectral photon radiance, [p/sec-cm 2 -sr-jum] obtained 

from Equation (3-7); 

T bg = the background temperature [K] of any radiant energy source other than 

the ground resolution element; and 
7 " G p = optical transmission [nd] of the cold filter. 

The effective solid angle fi (Figure 5) through which the detector receives energy from the ground 
resolution element is given by 

a = [sr] . - (5-27) 

4f 2 N 

The background effective solid angle f2 BG includes the effective solid angle £2 plus a little more for 
tolerance purposes. Ideally, they would be identical. 
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Solid Angles Effective Solid Angles 

0' = 27t(1-cos0) S2 = rain 2 ( 6 ) = — — 

0'bg= 27t(1-cos 6 d(?) 4f2 
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Figure 5. Background Viewing Geometry 





The effective solid angle of the background fi BG is given by 

£2 B g = *F S t sr ]> (5-28) 

where the view factor F<- is defined as 

Fg = F c for a cold shielded detector viewing the scene through a circular aperture 
(Appendix D), and 

F s = F a for a detector array surrounded by a cold fence (Appendix E). 


It follows that: 


F 


C 


= sin 



[sr] , 


(5-29) 


where 0 BG is the full-cone angle [deg] of the background. 


From Equation (5-20) one can write 
[A D Af]‘/2 
°* s NEP 


[cm-Hz^/W] 


and 


NEP 


R(X) 


vt 


Np + Nj„ 


1/2 


[W] 


Substituting Equation (5-31) into Equation (5-30) and setting N QD = 0 yields 

R(X)t, (A d Af) 'I 2 


D *BLIP( X ) “ 


qN P 


[cm-Hz 1/2 /W] 


(5-30) 


(5-31) 


(5-32) 


where D* B lip is the background-limited photon (BLIP) value of D*. 


Substituting Equations (2-3), (5-21), (5-23), and (5-24) into Equation (5-32) gives 


d *blip ( x ) 


X 

h c 



* ^2 

E'(X) dX + 

Jh 



- 1/2 

[cm-Hz^ 2 /W] . (5-33) 


Equation (5-33) includes only the photon noise and is for a photovoltaic detector. For a photo- 
conductive detector there is an additional term due to generation-recombination noise, and it re- 
duces D* BLI p by a factor of 
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5.2.2 Infrared System Noise 


5.2.2. 1 Quantization Noise 

The quantization noise is given by (Montgomery, 1978, p. B-l) 


S' 


N Q = 


SAT 


[e], 


(12) 1/2 2 q 

where Q is the number of bits used in the A/D converter. 
The infrared saturation signal is given by 

! A 2 

E'sat(^) dA [e], 

A, 


(5-34) 


(5-35) 


where the infrared saturation spectral photon irradiance E' SAT (\) is given by (see Appendix A) 


T o* 


E'sat(^) = L'sat(^) [p/sec-cm 2 -/zm] 


4f2 


(5-36) 


and where the infrared saturation radiance L , ^ A j(X) is the appropriate value to just saturate the A/D 
converter. 

5.2.2. 2 Other Infrared System Noise 

When the system noises are from unknown sources, they are designated as other system noise. 
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6 . FIGURES OF MERIT 


One important figure of merit that we will discuss here is the signal-to-noise ratio SNR, which is 
easily obtained by using the results of Section 2 for the signal S' and Section 5 for the noise N T0T . 
We will also discuss a figure of merit used in the visible and SWIR bands, the Noise Equivalent Delta 
Reflectance (NEAp), and for the infrared we shall discuss the Noise Equivalent Delta Temperature 
(NEAT). 

6.1 Noise Equivalent Delta Reflectance 

When visible sensors are assessed with respect to surface observations such as reflectance, the SNR is 
not always a convenient figure of merit. Users of space or airborne remote sensor data are frequently 
concerned with the characterization of ground targets through the measurement of variations in target 
reflectance. Because the variations of interest are often small in magnitude and are difficult to mea- 
sure precisely, there is considerable interest in the definition and measurement of the capability of the 
sensor to respond to small reflectance changes. This capability, related to sensitivity, is often described 
in terms of Noise Equivalent Delta Reflectance (NEAp), which is the minimum detectable variation 
in reflectance, and is sometimes preferred by the science user community over the spatial resolution of 
the system. 


For the visible and SWIR bands, NEAp is the amount by which p would need to change to cause the 
signal to change by an amount equal to the noise, or it is the smallest change in reflectance between 
two adjacent surface elements that can be resolved by the sensor. In this section we compute NEAp 
from the SNR. 

The figure of merit NEAp is given by 


NEAp = 


(I) S) 


s 

N 


where 7 (see Appendix F) is 
7 n 


7 = 


r AN 


(sec (j> - 1) 


[nd] 


[nd]. 


( 6 - 1 ) 


( 6 - 2 ) 
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and where 


To 



[nd] . 


(6-3) 


Also, <[>' is the angle between the line of sight and the surface normal (see Figure 6) and L A and L N 
are the atmospheric and scene radiances respectively, from the tables in Appendix B. 


The atmospheric optical transmission in the nadir direction r AN is given by 


r AN 



[nd], 


(6-4) 


where the total nadir optical thickness is obtained from the tables in Appendix B. 


6.2 Noise Equivalent Delta Temperature 
For the infrared bands NEAT is given by 


where L' 

t A 

s 

N 

dL's/dTs 



(6-5) 


= the scene photon radiance [p/sec-cm 2 -sr] , 

= atmospheric transmission [nd], 

= signal [e], 

= noise [e], and 

= the differential change in surface radiance with respect to surface temperature 
[p/sec-cm 2 -sr-K] . 


(See Appendix G for a detailed derivation.) 
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7. MODULATION TRANSFER FUNCTION (MTF) 


If the sensor were to scan a very low spatial frequency (sinusoidal variation in radiance) an amplitude 
variation AS(0) in signal would result. At spatial frequency k [cycles/mni] an amplitude variation of 
AS(k) would be obtained. The ratio of AS(k) to AS(0) is the Modulation Transfer Function (MTF) 
of the sensor; i.e., 

AS(k) 

MTF = nd (7-1) 

AS(0) 

7.1 Total MTF 

For a linear system, the total modulation transfer function is the product of the modulation transfer 
functions of the individual elements of the system. There are many different MTFs associated with a 
sensor system, which typically might include: 



mtf oa = 

Optical Aperture MTF 


mtf da = 

Detector Aperture MTF 


mtf sm = 

Satellite Motion MTF 


mtf cd = 

Charge Diffusion MTF 


mtf ct = 

Charge Transfer MTF 


mtf sj = 

Satellite Jitter MTF. 


The total or system MTF is the product of the component MTFs; i.e., 
MTF = MTF qa • MTF da • MTF sm • MTF CD • MTF CT 


7.2 Component MTFs 

In this section we list the MTF equations and their reference sources. 
7.2.1 Optical Aperture MTF 

The diffraction MTF is given by (O’Neill, 1955 and 1956) 


MTF 


OA 


(1 - 0 2 ) 


where 


( 7 - 2 ) 
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and where 


C 


m 


= -2j3 2 


for 0 < — < 
9 


Cl - £) 
2 


(7-7) 


2 . 2/3 . / 1 


C m = -2/3 2 + Sini// + 

111 7T 


* 


-2 I - — \ Tan" 1 


1 + £ 

l - p 


Tan 


* 


1 - j3 

for — < 


A 1 + /3 
< (7-8) 


and 


C 


m 


= 0 



1 + P 


(7-9) 


and \p is given by 


= Cos 


-1 


1 + /3 2 - A 2 
2/3 


with 



where k is the spatial frequency [cycles/mm] measured in the image plane and 


(7-10) 


(7-11) 
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k 0 = ~^T~ [cycles/mm] . 
2At N 


(7-12) 


The modulation transfer function MTFq^ = 0 for A = 2 at the cutoff frequency when k = k^, where 
k^ = 2 kg [cycles/mm], (7-13) 


and where 


f N = — [nd] . 


(7-14) 


where X is the wavelength, f is the focal length of the optics, and D is the diameter of the optics; and 


0 = — [nd] 


(7-15) 


where D 0 is the diameter of any obscuration. 


Note that the quantity k Q given in Equation (7-12) is not defined in O’Neill’s paper but must be de- 
fined this way to be consistent with his Figure 3 in which the MTF goes to zero at A = 2.0. 


7.2.2 Detector Aperture MTF 


The detector aperture MTF is given by (Jensen, 1968, p. 27) 


sin(7rkd s ) 


(7-16) 


where dc 


the detector width [mm], and 

spatial frequency in the image plane [cycles/mm]. 


7.2.3 Satellite Motion MTF 

The MTF due to linear image motion is given by (Jensen, 1968, p. 

sin(7rkVjtj) 

MTFc M = 

SM TTkVjtj 


(7-17) 


where the image velocity Vj is given by 


[ km/sec ] , 
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where V SUB = the subsatellite point velocity [km/sec] ; 

k = spatial frequency in the image plane [cycles/mm] ; and 

tj = integration time [sec]. 


7.2.4 Charge Diffusion MTF 

The charge diffusion MTF is given by (Jespers, 1975, p. 519) 

exp - a a d 


mtf cd = 


1 + <*a L 
exp - a a d 


1 + a a L o 

where d = photodetector depletion region depth (typical value = 5 jam). 


(7-19) 


The silicon absorption coefficient a is a function of wavelength and temperature and is given by 


a = 10 z [cm -1 ] 


(7-20) 


where, after curve fitting to Jespers’ Figure 25 for silicon, we have 


z = 2.897652 - 4.044143 (X - 0.82) - 5.219219 (X - 0.82) 2 
- 3.828495 (X - 0.82) 3 + 22.16724 (X - 0.82) 4 , 


(7-21) 


where X = 


L 0 = 

and where 


wavelength [jum], 

diffusion length (typical value = 50 /am), 


L = 


L 2 
L 0 


1 + (27rkL 0 ) 2 


1/2 


(7-22) 


7.2.5 Charge Transfer MTF 

The MTF due to inefficiency in charge transfer in a CCD detector is given by (Jespers, 1976, p. 520) 


MTFfT e 


-M CT e 


1 - cos 


7rk 


max 


(7-23) 
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where e is the charge transfer inefficiency and where the number of charge transfers Mqj is given by 
M C t = m s m p , (7-24) 

where m s = the number of stages, detectors, or picture elements, and 

m p = the number of clock phases for readout. 

The Nyquist spatial frequency k max is given by 

^ma x rp 5 (7-25) 

where P D is the detector pitch. 

7.2.6 Satellite Jitter MTF 

The satellite jitter MTF is given by (Jensen, 1968, p. 124) 


MTFsj = J 0 (27rkf© M ) 


where 0 


M 


maximum satellite angular movement [rad], 
zeroth order Bessel function, and 
focal length (mm). 


(7-26) 
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APPENDIX A 


DETECTOR IRRADIANCE 


The objective of this appendix is to show the relationship between the scene spectral radiance and the 
irradiance on the detector. 


When a sensor images an area on the surface of the Earth called the instantaneous field of 
view (IFOV), it also receives energy from the intervening atmosphere. The scene spectral radiance 
L(X) [W/cm 2 -sr-jum] is defined as the combined spectral radiance from the atmosphere and the 
IFOV area A, as viewed from the sensor. The power that the area Aj and the intervening atmosphere 
radiate through the solid angle into the sensor entrance aperture and into the detector is then 


f *2 

4) = d0(A) [W] 

h i 


(A-l) 


where 


d0(A) = r 0 FZ 0 A[L(A) dA [W] (A 2) 

and where t q = the sensor optical transmission. However, by definition, for a detector of 
area A D , Equation (A-2) may also be written in terms of the detector spectral irradiance E(A) 
[W/cm 2 -pm] as 

d 0(A) = A D E(A) dA [W] (A-3) 

By comparing Equations (A-2) and (A-3), one can express the scene spectral irradiance E(A) in terms 
of the scene spectral radiance L(A) and sensor design-related parameters A D , Aj, and as 

t £2 Ai 

E(A) = L(A) [W/cm 2 -pm) . (A^) 

A D 


Now, from Figure A-l, 




D 

f 

H 


a 


0 

ft 

ft 

ft 

ft 


0 

1 

d 


= Detector area (/urn 2 ) 

= Area of sensor entrance aperture (cm 2 ) 

= Area of ground resolution element (cm 2 ) 

= Diameter of sensor entrance aperture (cm) 

= Effective focal length (cm) 

= Satellite height (km) 

= Instantaneous angular FOV (as in previous figs) (rad) 

= Half-cone angle of optics (deg) 

= The effective solid angle through which the detector receives energy (sr) 

= Solid angle subtended by the sensor entrance aperture at the ground resolution element (sr) 

= Solid angle subtended by the ground resolution element at the sensor entrance aperture (sr) 

= Solid angle subtended by the detector at the entrance aperture (sr) 


Figure A-l . Radiometric Geometry 
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where a 


the sensor instantaneous field of view (IFOV), and 


H = the distance from the satellite to the area Aj. 


Also, for a square detector 


a d = « 2f2 > 


(A-6) 


where f is the sensor optics effective focal length, which represents the optical path length and is ap- 
proximately equal to the geometric focal length for small optical convergence angles (<C10 ), and 


ttD 2 


a = 

0 -j2 


H z 4H Z 
where D = aperture diameter of the sensor. 


(A-7) 


Using Equations (A-5), (A-6) and (A-7), one can write 


A i n o n 


V D 4f d 


N 


where the sensor optics f-number is given by 


f N " “d ' 


The half-cone angle d of the optics is given by 

D 


sin 6 = 


It is also given by 


tan 6 


2f 


D 

IF 


The f-number f N is related to the approximate f-number 


f N - 


f ' , 

) by the expression 


f' 

d" 



2 I 


1/2 


Substituting Equation (A-8) into Equation (A-4) (See Figure A-2) gives 


(A-8) 


(A-9) 


(A-10) 


(A-ll) 


(A- 12) 


A-3 




A(%) 

1 

13.4 

1.1 

10.9 

1.2 

9.1 

1.3 

7.7 

1.4 

6.6 

1.5 

5.7 

2 

3.2 

3 

1.9 

4 

0.8 

5 

0.5 


f N = F-NUMBER 
f 1 



A =(V^-)/ 100, N 


Comparison of f N 
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Figure A-2. 


and F /D 



E(X) = 


L(X) [W/cm 2 -pm] . 


(A- 13) 


7TT„ 


4f 


N 
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APPENDIX B 


COMPUTED EARTH ATMOSPHERE RADIANCES 


With the permission of the authors, the following discussion has been taken from a NASA unpublished 
report (Mattoo, 1984) and a memorandum (Fraser, 1981). 

The computed radiances at the top of the Earth-atmosphere system are tabulated for 15 wavelengths. 
The computations are made for plane-parallel models of the Earth-atmosphere system; that is, the op- 
tical properties of the models do not vary in a horizontal plane. However, the vertical profiles of the 
concentrations of the atmospheric constituents are arbitrary, but realistic. The ground reflects light 
according to Lambert’s law, which implies that the radiance of the reflected light is constant, indepen- 
dent of direction. 


Only one atmospheric model is used, and it contains the standard dry gas and the gases that absorb in 
the spectral bands of interest. The variable trace gases are 316 Dobson units of 0 3 and 2.5 cm of 
H 2 0. The model also contains particulates, but not liquid or ice clouds. 

The average normal optical thickness of each constituent is given for each spectral band in Table B-l. 
The total nadir optical thickness (5 TN ) at a wavelength equals the sum of the optical thicknesses of the 
constituents: 


5 TN " 5 R + 5 G + 5 A ’ 

where 5 R = the scattering optical thickness (Rayleigh) of the dry atmosphere, 

5 a = the optical thickness of the aerosols (particulates), and 

5 g = the optical thickness of the absorbing gases. 

The total optical transmission r for the direct sunlight is given by 
E Xs 

t - = exp(-5 TN sec 8 Z ) 

E Xo 


(B-l) 


(B-2) 


where E^ s and E^ Q are the spectral irradiance of the direct sunlight at sea level and above the atmos- 


B-l 


Table B-l. Optical Thicknesses for Scattering by Molecules (6 R ), Absorption by Gases (5 G ), 
and Scattering from Aerosols (S A ) 


X 

0.4 

0.44 

0.48 

0.52 

0.56 

5 r 

0.3637 

0.2451 

0.1713 

0.1234 

0.0912 

5 g 

0.0000 

0.0007 

0.0048 

0.0151 

0.0296 

5 a 

0.3600 

0.3273 

0.3000 

0.2796 

0.2571 

5 tn 

0.7237 

0.5731 

0.4761 

0.4181 

0.3779 

X 

0.62 

0.66 

0.70 

0.74 

0.82 

5 r 

0.0603 

0.0468 

0.0368 

0.0294 

0.0195 

5 g 

0.0325 

0.0169 

0.0135 

0.0145 

0.0560 

5 a 

0.2370 

0.2228 

0.2087 

0.1988 

0.1794 

5 IN 

0.3298 

0.2865 

0.2590 

0.2427 

0.2549 

X 

0.88 

1.05 

1.25 

1.60 

2.20 

5 R 

0.0147 

0.0072 

0.0036 

0.0013 

0.0004 

5 g 

0.0039 

0.0000 

0.0092 

0.0066 

0.0608 

5 a 

0.1694 

0.1420 

0.1193 

0.0932 

0.0682 

5 tn 

0.1880 

0.1492 

0.1321 

0.1011 

0.1294 
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phere, respectively, and Q z is the solar zenith angle. The values of E^ 0 are given in Table B-2. The 
quantity E^ 0 is the solar irradiance on a square centimeter of surface perpendicular to the solar rays. 

The aerosols (particulates) are assigned properties of those occurring in continential regions. The size 
distribution function (n) of the particle radius (r), which is the number of particles per cubic centi- 
meter of air per micrometer of radius, decreases very rapidly with increasing radius: 


The index of refraction of the particulates is m = 1.4300 - 0.003 5i. The nadir radiances of the model 
are not sensitive to the vertical profile of the particulate concentration. Here, a realistic profile with a 
high concentration near the ground is assumed. 

The computations are made with a computer code developed by Dr. J. V. Dave of IBM as modified 
by R. S. Fraser of NASA. The equation of radiative transfer is solved numerically by a procedure 
that iteratively accounts for successive scatterings of light from the atmosphere and the ground. I he 
polarization characteristics are not accounted for, and as a result, the computed radiances are in error 
by a few percent. 


The input parameters for the computations include the model parameters: the vertical profiles of the 

concentrations of the gases and aerosols (particulates), the scattering and absorption optical thickness, 
the gaseous absorption coefficients, and 10 values of the surface reflectance p. The volume extinction, 
scattering and absorption coefficients and scattering phase function of the particulates are computed 
according to the Mie theory by a separate code, and are part of the input. 


The nadir spectral radiance at the top of the atmosphere L N can be expressed as follows: 

-5, 


TN 


N 


L N = L s e + L a , 


(B-4) 


where the first term on the right-hand side of the equation gives the radiance at the surface (ground) 
L s (assumed to be Lambertian) attenuated by the atmosphere; the second term gives the radiance of 
just the atmosphere L^, or path radiance. The radiances L N and are given in Table B-2, where 


B-3 



>2 


X = the center wavelength [mm], 

6 Z = the solar zenith angle [rad] , 

p = the surface reflectance [nd], 

5j N = the total optical thickness [nd], 

E^ q = the solar spectral irradiance [mW/cm 2 -|um] , 

L n = the total nadir spectral radiance [mW/cm 2 -sr-/am] and 

= the atmospheric path spectral radiance [mW/cm 2 -sr-pim] . 
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Table B-2. Reflected Solar Spectral Radiance (Continued) 
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Table B-2. Reflected Solar Spectral Radiance (Continued) 
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Table B-2. Reflected Solar Spectral Radiance (Continued) 
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Table B-2. Reflected Solar Spectral Radiance (Continued) 
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Table B-2. Reflected Solar Spectral Radiance (Continued) 
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Table B-2. Reflected Solar Spectral Radiance (Continued) 
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APPENDIX C 


SENSOR FIELD OF VIEW 


Each time the satellite completes an orbit around the Earth, it maps out a swath on the Earth’s sur- 
face of width S w . The angle subtended at the satellite by S w is the sensor field of view O. Referring 
to Figure 1 and using the law of cosines, one can write 


0 = 2 Cos -1 



+ S 2 d - R 2 

2S d R s 



[rad] , 


(C-l) 


where 


Rg = H + R e [km] 


(C-2) 


and 


S d = 


1 1/2 


R 2 + R 2 - 2R s R e Cos f-y 


[km] 


(C-3) 


From Figure 1, it can also be seen that 


<Pc = 


~rT 


[rad] . 


(C-4) 


The swath width S^ is given by 


S w - 


Ad c 

“nT 


[km] 


(C-5) 


where d(- is the cross-track distance covered by the sensor and, for total coverage at the equator, is 
given by 

d c = 27rR e [km] . (C-6) 


Now, if the total time required to map the Earth t MA p is known, then the total number of swaths Ng 
required to map the Earth is 


N c = 


l MAP 


[nd] , 


(C-7) 


where tg is the satellite’s orbital period or the time required to map one swath and is given by 


C-l 



where S' Q is the percentage overlap across the ground track. 


C-2 


APPENDIX D 


VIEW FACTOR FOR A SINGLE DETECTOR VIEWING A CIRCULAR BACKGROUND 

The view factor F c for a detector that is receiving radiation from a circular background can be de- 
fined as 


F 


C 


^d^aa 


[nd] , 


CD-I) 


where <F' = the photon flux into the detector [p/sec] , 

A d = the area of the detector [cm 2 ], 

and 

L AA = I B'(X,T B g) dA [p/sec-cm 2 -sr] , (D-2) 

* 0 


where B'(A,T bg ) = Planck’s function evaluated at the background temperature T bg , and 
X c = the detector cutoff wavelength. 

If a detector of area A D views a source of radiance L'^ with area dA' through a solid angle ft, the 
flux into the detector is 

d<F' = L'^ft Cos(0' c ) dA' [p/sec] . (D-3) 

The angle 4>' c is shown in Figure D-l. 

However, 

dA' = 27rr' dr' [cm 2 ] (D-4) 


and 

A d Cos (</>') 

ft = — [sr] . 

S 2 

Substituting Equations (D-4) and (D-5) into Equation (D-3) gives 


(D-5) 


D-l 



a = Circular background radius (cm) 

A d = Area of detector (jum) 

d = Distance between the detector and the circular background (cm) 
dA' = 27rr'dr' = Area of elemental annular ring of radius r' (cm 2 ) 

S = Distance from center of detector to elemental area on circular background (cm) 

0 c = Half angle subtended by the circular background at the center of the detector (deg) 

<p' c = Half angle subtended by the elemental circular area dA' at the center of the detector (deg) 

dn = Elemental solid angle subtended by the detector at a point on the elemental area dA' (sr) 


Figure D-l. View Factor Geometry for a Single Detector Viewing a Circular Background 


D-2 



d<F' = 


27rrA D L' AX Cos 2 (0' c )dr' 


(D-6) 


However, from Figure D-l we see that 

d d 


C0S = S ( d 2 + r '2 ) 1/2 


[p/sec] . 


[rad] , 


and substituting Equation (D-7) into Equation (D-6) gives 


d<F' = 2?rd 2 A d L' ax 


r' dr' 


(d 2 + r' 2 ) 


[p/sec] . 


It follows that 


4>' = 27rd 2 A 


! a 

r ' dr ' 

„(a 2 ♦ r ' 2 ) 


[p/sec] . 


However, 


r' dr' 


(d 2 +r' 2 ) 2 2d 2 2(d 2 + a 2 ) 


Therefore, substituting Equation (D-10) into Equation (D-9) gives 


$' = ttA d L' ax I 1 - 


d 2 +a 2 


[p/sec] 


However, 


(d 2 + a 2 ) l,! s "" c 


= = Cos <p r [nd] 

'a 


Therefore, by substituting Equation (D-l 2) into Equation (D-ll) we obtain 


$' = ttA d L' ax Sin 2 <f) c [p/sec] 


(D-7) 


(D-8) 


(D-9) 


(D-10) 


(D-ll) 


(D-l 2) 


(D-l 3) 


D-3 



APPENDIX E 


VIEW FACTOR FOR A DETECTOR IN AN n X m ARRAY VIEWING 
A RECTANGULAR BACKGROUND 


The view factor F A for a detector in an n by m array viewing a rectangular background (Figure E-l) 
is given by 


F 


A 


<fr' 

^D^AA 


(E-l) 


where <F’ is the total number of photons per second entering the detector from the background through 
a rectangular aperture, given by 


<F' — \ Cos dAg [p/sec] 

'0 


and, as in Appendix D, is 

L'ax = \ B’(X,T bc )dA [p/sec-cm 2 -sr] 

Jo 


(E-2) 


where the photon radiance is given by Equation (E-2) and 


where B'(A,Tg G ) 


£2 d 

dAg 

4>f 


the Planck function evaluated at the background temperature T BG 
[p/sec-cm 2 -sr-pm] ; 

the solid angle subtended by the detector at an arbitrary point on the rec- 
tangular aperture through which the background is viewed; 
the differential area on the background; and 

the angle between the normal to the differential area dAg and the line be- 
tween the center of the detector and the center of the differential area dAg. 


If D is the distance between the center of the detector and the differential area dAg then 


E-l 



A b = Background area (cm 2 ) 

A d = Detector area (firvi 2 ) 

D = Distance from center of detector to dA b (cm) 
dA b = Differental background area (cm 2 ) 

H f = Fence height (cm) 

Lf = Fence length (cm) 
p = Detector pitch (cm) 

W f = Fence width (cm) 
d = Optics half-cone angle (deg) 

(f> f = Differential area normal angle 

Figure E-l. Geometry for a Detector in n X m Array Viewing a Rectangular Background 


E-2 



A d Cos($ f ) 


[sr] , 


(E-3) 


i2 d _ 


D J 


where 


Cos%) = [nd] . (E-4) 

To evaluate the quantities D, <j> f , and z in Equations E-3 and E^l, it is convenient to define two sets 
of coordinates. The two systems (x,y,z) and (x',y',z') are shown in Figure E-l. The n X m detector 
array lies in the x-y plane and the plane of the rectangular aperture, through which the background is 
viewed, is a distance Hj- from the x-y plane. 


Consider an n X m array of detectors, and let p be the pitch (distance between detector centers) in 
both directions (Figure E-l). Let the detector array be symmetrically surrounded by very cold sides 
(fence) of length L f , width W f , and height H f . The length and width are adjusted for a given height 
so that the edge detectors in the array can just accommodate the optical bundle which is defined by 
half the bundle cone angle 6. The relationships for the length and width W f of the fenced area 
that surrounds the detector array are 


L f = np + 2H f Tan(0) [m] , 


(E-5) 


and 

W f = mp + 2H f Tan(0) [m] , (E-6) 

where n = the number of detectors along the length of the detector array, and 

m = the number of detectors along the width of the detector array. 


Equations (E-5) and (E-6) follow because 


Sin(0) = • 

/T n 

The coordinates of detector i, j with respect to the edges of the fenced area are given by 


(E-7) 


E-3 


X d = H f Tan(0) - — + i 


IP 


(E-8) 


and 


where i = 

j = 


y d = H f Tan(0) — + jp , 

the number of detectors along the fence length, (i = 1, n) and 
the number of detectors along the fence width. O’ = 1, •••, m). 


(E-9) 


The coordinates of the arbitrarily placed differential area dA B on the surface formed by the top edges 
of the fence are (x,y,z). The coordinates with respect to the center of the detector i, j are given by 


and 


where 


and 


X = X - X d , 


y = y - y d 


z = z , 


z = H r 


D = (x' 2 + y' 2 + z' 2 ) 


(E-10) 


(E-ll) 


(E-l 2) 


(E-13) 


(E-l 4) 


By substituting Equations (E-10) through (E-13) into Equation (E-14), one gets 


D = 


( x - x d ) 2 + (y - yd ) 2 + z2 


Letting 


dA b = dx dy 


(E-l 5) 


(E-l 6) 


and substituting Equations (E-3), (E-4), (E-13), and (E-l 6) into Equation (E 2) results in 

dx dy 


*’ - a d H f 2 B' ax 



[ p/sec ] . 


(E-l 7) 





Substituting Equation (E-17) into Equation (E-l) yields 


F 


A 



71 


Wf 

^ g(y) dy , 
0 


where 



Equation (E-l 9) may also be written as 



where 


(E-l 8) 


(E-l 9) 


(E-20) 


X = ax 2 + bx + c 

(E-2 1 ) 


a = 1 

(E-22) 

7- 

b = -2x d 

(E-23) 

= 

and 



c = x^j + (y - y d ) 2 + z 2 . 

(E-24) 


It follows from Equations (E-22) through (E-24) that 



4ac - b 2 = 4(y - y d ) 2 + 4z 2 . 

(E-25) 

= 

Now, Equation (E-25) reveals that 


r^T 

4ac > b 2 . 

(E-26) 



Therefore, Equation (E-20) may be written (Dwight, 1947, p. 33, Equations 160.01 and 160.02) as 


E-5 


II 11 I 



4a 


2ax + b 


(E-27) 


dx 2ax + b 


f dx _ 

lit " 


Tan 


rX 


Ro 


R, 


where 


and 


and 


r = 4ac - b 2 


Rj = r = ( 4ac - b 2 


R 2 = rRj = r 3/2 • 


1/2 


(E-28) 


(E-29) 


(E-30) 


Applying Equation (E-27) to Equation (E-28), one gets 


g(y) 


2aL f + b 


r ^aL 2 + bL ( - + ^2 


4a , l 2aL f + b 

+ Tan -1 


R 


b 4a 

r c Rt 


Tan' 


R, 


(E-3 1 ) 


Use of Equation (E-31) in Equation (E-18) and numerical integration enable the form factor F A to 
be computed. 


E-6 


APPENDIX F 


COMPUTATION OF 7 


Substituting Equations (2-11) and (2-12) into Equation (2-9) yields 


S' = tjT, 


■**(■*) fe) J “ 


L(X) dA [e] 


Therefore, the signal equation for the visible and SWIR bands is given by 


(F-l) 


S ' = VqAd 7 ? 




LAA [el 


(F-2) 


where we have replaced the integral in Equation (F-l) by LAA because L(A) varies slowly over the 
spectral bandpass AA. Taking the differential of Equation (F-2) with respect to the scene radiance, 
one obtains 


dS ' “ ^oAd 77 


£) fe) 


dL AA [e] 


(F-3) 


Dividing Equation (F-2) by (F-3) gives 


S 

n" 


where the signal S = 


S' 

dS 7 

S', 


L 

dL 


[nd] 


and the noise N = dS 


(F-4) 


However, 



[W/cm 2 -sr-pm] . 


Substituting Equation (F-5) into Equation (F-4) gives 


S 

N 



[nd] . 


(F-5) 


(F-6) 


F-l 


Letting 


7 = 


/ dL 

\dp 


[nd] , 


dp = NEAp 

and substituting Equations (F-7) and (F-8) into Equation (F-5), one obtains 


(F-7) 


(F-8) 


However, 


N NEAp 


L = L§ + L» [W/cm 2 -sr-pm] 


where the surface spectral radiance Lo has the functional form 


L| = pK| [W/cm 2 -sr-pm] 


and the atmospheric spectral radiance is a constant with respect to p and is given by 


La = ^2 [W/cm 2 -sr-pm] 


(F-9) 


(F-10) 


(F-ll) 


(F-12) 


The superscript S denotes that the sensor observes these spectral radiances along a slanted path. (See Fig- 
ure 6.) Since the scene spectral radiances are assumed to be equal in the normal and slant directions, no 
superscript is used. 

Substituting Equations (F-l 1) and (F-12) into Equation (F-10) gives 


L = KjP + K 2 [W/cm 2 -sr-pm] . 


Differentiating Equation (F-l 3) with respect to p results in 


= K, [W/cm 2 -sr-jum] 

dp 1 


Solving for Kj in Equation (F-l 1) and substituting it into Equation (F-14) results in 

T s 

dL 


= — [W/cm -sr-pm] . 
dp p 


(F-l 3) 


(F-14) 


(F-l 5) 


F-2 


Hill I I II; I 


Substituting Equation (F-15) into Equation (F 7) gives 


7 = 




p [nd] . 


(F-16) 


The spectral radiances L and L| are observed along the line-of-sight direction with angle of 0 and a 
line-of-sight surface-normal angle <1>’ (Figure 6). These angles are related by 

-1 


0' = Sin 


[K 


sin <p 


[rad] . 


(F-l 7) 


The surface spectral radiance Lg is given by 


Li = 


7T 


pr AN^ ^ [W/cm 2 -sr-Mm] 


(F-l 8) 


where E s = the irradiance at the surface of the Earth [W/cm 2 -pm] 

The atmospheric transmission along the nadir direction is given by 
t an = e 5 ° [nd] , 

where 6 0 is the optical depth along the nadir direction. (See Appendix B) 

The surface spectral radiance L s observed along the nadir direction (for which <p - 0) is given by 
E S 

Lg = pr AN [W/cm 2 -sr-pm] . 

Dividing Equation (F-l 8) by Equation (F-l 9) gives 

L S = L S T i S N C 4)1 “ n [W/cm 2 -sr-pm] . 

Substituting Equation (F-21) into Equation (F-16) we obtain 


(F-l 9) 


(F-20) 


(F-21) 


- 7 o 

r = d s N c *' - 0 


[nd] 


(F-22) 


F-3 


where 


To = 


Lp 

T N 


[nd] . 


(F-23) 


We assume the total spectral radiance L is the same along the nadir and slant directions; hence 


L = Lg + [W/cm^-sr-pm] 


(F-24) 


or 


Lo = L - La [W/cm 2 -sr-pm] . 


(F-25) 


Substituting Equation (F-25) into Equation (F-23) gives 


T 0 = 


L - L 


N 

A J 


p [nd] 


(F-26) 


or 


T 0 - F 


N 


[nd] 


(F-27) 


1 - 


F-4 


APPENDIX G 


NOISE EQUIVALENT DELTA TEMPERATURE 


Equation (2-9) enables the signal to be written as 
S' = tjAprjE'^ [e] , 

where 

C * 2 

E'^ = \ E'(X) dX [p/sec-cm 2 ] . 

X 1 


However, 


7T7 0 

E'(X) = L'(X) [p/sec-cm 2 -jum] 


4P 


N 


and 


L'(X) = t a B'(X,T s ) + e A B'(X,T A ) [ p/sec-cm 2 -sr-//m] . 


Therefore, 


E'ax * — — L' ax [p/sec-cm 2 ] 

« 2 n 

where 

L'^ a = 1 L'(X) dX [ p/sec-cm 2 -sr] . 

X 1 

Equation (G^I) can be written as 


L AX = r A L 'sAX + e A L AAX [p/sec-cm 2 -sr] , 

where X 2 

L'sax = ^ B'(X,T s ) dX [p/sec-cm 2 -sr] 

X 1 


G-l 


(G-l) 


(G-2) 


(G-3) 


(G-4) 


(G-5) 


(G-6) 


(G-7) 


(G-8) 


and 


l 'aAX “ 



B'(X,T a ) dX 


[ p/sec-cm 2 -sr] . 


Taking the differential of Equation (G-l), one obtains 

dS = (tjAj-jTj) dE ^ [e] . 

Replacing dS' by N, and dE'^ A by NEI results in 

N = (tjAjjP) NEI [e] , 

where N = the total noise [e], and 

NEI = the noise equivalent photon irradiance into the detector [p/sec-cm 2 ]. 

Taking the differential of Equation (G-5) gives 

7TT 0 

dE AX = dL AX [p/sec-cm 2 ] . 

4f 2 N 


(G-9) 


(G-10) 


(G-l 1) 


(G-l 2) 


Replacing dE'^ A by NEI, and dL'^ A by the noise equivalent photon radiance (NEPR) gives 


7 TT„ 


NEI = NEPR [p/sec-cm z ] . 


(G-l 3) 


4P 


N 


Taking the differential of Equation (G-7) gives 

dL 'AX = r A dL SAX + 6 A dL A AX [ p/sec-cm 2 -sr] . (G-14) 

However, since we are not interested in perturbations due to changes in the atmosphere, we assume 
L'aAX is to be constant - Equation (G-14) then becomes 

cIL'ax = T a dL SAX [p/sec-cm 2 -sr] . (G-l 5) 


G-2 


Replacing dL' AX with NEPR, and dL' SAX with the surface NEPR, NEPR S in Equation (G-15), gives 


NEPR = r A NEPR s [p/sec-cm 2 -sr] . 


(G-16) 


However, by definition, 


or 


dL' 


dL' 


SAX 


SAX 


dTc 


dTc 


[p/sec-cm 2 -sr] 



Replacing dTg with NEDT and dL^ AA with NEPRg in Equation (G-18) yields 


(G-l 7) 


(G-18) 


and 


NEAT = 


dL 'sAX 

dT s 



[p/sec-cm 2 -K] , 


where the Planck function is given by 


(G-l 9) 


(G-20) 


B'(X,T S ) 



[ p/sec-cm 2 -sr-/rm ] . 


Differentiating Equation (G-21) with respect to T s yields 


dB'(X,T s ) 

dT s 


C 2 X 3 ^B'(X,T s )^exp 


C’.Ti 


[p/sec-cm 2 -sr-pm-K] 


(G-21) 


(G-22) 


G-3 


Equation (G-19) may be written in terms of signal-to-noise ratio (S/N) as follows. 
(G-l) by Equation (G-ll) gives 


s _ e 'aa 

N NET 


[nd] . 


Substituting Equations (G-5) and (G-l 3) into Equation (G-23) and replacing dL'^ 

s l 'aa 

— = nd . 

N NEPR 


Substituting Equation (G-l 6) into Equation (G-24) gives 


S L AX 


N 7 A NEPR S 


[nd] 


or 


L' 


NEPRg = 


AX 


(I) 


[p/sec-cm 2 -sr] 


and, finally, by substituting Equation (G-26) into Equation (G-19), one obtains 


L' 


AX 


' \/ dL 's 


Dividing Equation 

(G-23) 
with NEPR gives 
(G-24) 

(G-25) 

(G-26) 


NEAT = 


S 


[K] . 


(G-27) 
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